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Abstract 

We analytically calculate to leading order the asymptotic form of quasinormal frequencies of 
Kerr black holes in four, five and seven dimensions. All the relevant quantities can be explicitly 
expressed in terms of elliptical integrals. In four dimensions, we confirm the results obtained by 
Keshest and Hod by comparing the analytic results to the numerical ones. 
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1 Introduction 



Perturbation of black holes are known to reveal characteristic damped oscillation modes which domi- 
nate the time evolution in certain intermediate period of time [1]. Since the frequencies are complex, 
they are called quasinormal modes (QNMs). For a general review and classification, see Refs. [2, 3, 4]. 
They depend only on the fundamental parameters of the black holes, such as mass, angular momen- 
tum and charge. For a Schwarzschild black hole, an asymptotic formula for high overtones has been 
obtained using the monodromy matching method [5]: 



in the units G = c = K = 1. The analytic value In 3 in the real part of the above formula was used to 
argue that the relevant gauge group in loop quantum gravity should be SO (3) [6]. Although this result 
turns out to be not universal [5], one still expects QNMs to play an important role in understanding 
black holes and quantum gravity. 

The monodromy method has been generalized so that first order correction to the asymptotic 
form of quasinormal frequencies can also be calculated [7]. In the case of Schwarzschild black holes, 
comparisons with numerical results have been made and the agreement is excellent [8, 9]. In Ref. [10], 
the calculation is further extended to the second order. The agreement is not as good as the first order 
case. In fact, there seems to be sizable discrepancy in higher angular momentum case (1 = 6). More 
study is needed to clarify the situation. 

Although high overtones of QNMs in spherically symmetric black holes have been extensively 
studied, they are not as well studied in the case of rotating black holes. In Ref. [11], convergent 
numerical results have been presented. Making use of the monodromy analysis used in Refs. [5, 12], 
Keshet and Hod obtain the first analytic results that are in excellent agreement with the numerical 
ones [13] . There are two distinct properties of rotating black holes that makes its highly damped QNMs 
quite different from that of spherically symmetric ones. First, there are more than one turning points. 
Second, there is a term linear in the frequency uj, which contributes to the real part in leading order. 
In this paper, we will extend their method to find the asymptotic formula of quasinormal frequency 
for higher dimensional Kerr black holes. 



2 Calculation of the asymptotic form of quasinormal frequen- 



Consider a massless scaler field $ in the background of a D-dimensional Kerr black hole. Using the 
Boyer-Linquist coordinate, we can write the metric as 




(1) 



cies 




sm' 



2 



e 



ds 2 = -^ 



^ + de 2 + 



[adt - (r 2 + a 2 



)d<p] 2 + r 2 cos 2 9dfl 2 L 



D-4- 



(2) 



Here. 



p 2 = r 2 + a 2 cos 2 9, 
A = r 2 -pr 5 - D + 



(3) 
(4) 



2 



dCl 2 - ) _ 4 is the metric on SD-4/Z2, with Qd-4 = {i>i, ■ ■ ■ , ipD-i) an d < ipi < n for all i. The location 
of an event horizon is determined by A = 0. For D = 4, there are both the outer and inner horizons 
r± = /i/2 ± J /j, 2 / 4 — a 2 . For D > 5, there is only one horizon. The surface r — r + is also a Killing 
horizon. The corresponding Killing vector is £ = d t + Qh9<p- Here, = ^ ° a2 is the angular 

velocity on the horizon. From £, one can find surface gravity k = (D The ADM mass 



is given by M 



(D-2)A D - 2 y. 

167T ' 



with An 



27r'" +1 >/ 2 



the area of an n-dimensional unit sphere. The 



Bekenstein-Hawking entropy is given by S = \A D - 2 {r\ +a 2 )r^- 4 [14]. Making use of the identity 
dS = ^(dM - fl H dJ), one then find J = 



A D - 2 (ri+a 2 )a 



It's well known that <& satisfies the Klein-Gordon equation: 
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where u — cos 6 and Yz,,Mi,...,Mr>_5 the spherical harmonics. i?; TO z,(r) and Si m L{u) satisfy 
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respectively. A; m i is the separation constant. Define T/„,.l(u) 
angular equation becomes 
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When a — > 0, the solution to the above equation becomes the hypergeometric function 2 Fi(A, B: C, u) 
with 
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In such limit, Tj® L (u) together with e m * and Y L . Mi,...,m d _ 5 (Qd-4,) form the spherical harmonics on 
S D ~ 2 , and therefore A; m ^ = 1(1 + D — 3). As a result, A = |(m + L — /) = — p, where p specifies the 
number of zeros of T^ L (u) for u in the interval [0, 1]. 
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It is straightforward to generalize Flammer's method to arbitrary dimension D [15]. In the limit 

uo — ► — zoo, 



MrnL = A C+ Ai +0( C 1 ), 

Ti mL {u) = Z (x) + -Z x {x) + 0(c- 2 ), 
c 

where c = iuoa and a; = cu. To leading order in c _1 , Zo(a;) satisfies 

^'(s) + i(D - 3 + 2i)4(x) - l(a; - A )Z (x) = 0. 
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The solution is given by Zq(x) = c X / 2 L^. 



, , >L Ax), where L p n '(x) is the generalized Laguerre 
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polynomials with p zeros. For Z {x) to have the same number of zeros as T^ L (u), we find 

A = [2(1 - \m\) + (D - 3)]. 



(16) 



This is consistent with the result in the D = 4 case [11]. Let Ri m L(r) — [r D 4 A] 1 ^ 2 Ri m L(r), and the 
radial equation becomes 
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Similar to Ref.[13], we define z = fV(r')dr' with V(r) = A' 1 (r)[q (r) + u^q^r)} 1 / 2 . The radial 
equation becomes 
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It has been shown using the monodromy matching method that the condition for quasinormal modes 
is given by 
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Here, Ct,o is a contour in the complex plane that runs from t\ to ti lying beyond the outer horizon 
r = r+. ti,*2 are the transition points in the right half complex plane determined by the condition: 



q o {t)+uj- 1 q 1 {t)=0. 



(24) 
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Figure 1: The Stokes lines (dashed) and anti-Stokes lines (solid) emanting from the turning points t\ 
and ti in complex r-plane for D = 4 a = 0.5 are shown. r_ and r+ are the inner and outer horizon 
radii, repectively. 

In the limit lu — > — ioo, 



2c/ Vdr ^2u> I ^±dr+[ -^Ldr. (25) 
Making use of the explicit expression for qi(r) and following the convention in Ref. [13], we have 



So = -2* I ^dr; (26) 

5 m = 2i f M ° r !JL dr; (27) 

S A = 2i [ — % dr. (28) 

The sign convention is opposite to that in Ref. [13]. It is chosen to make 5a,5 m and S\ all positive. 
The QNM quantization condition in eq (23) becomes 

luSq + m6 m + iAo5\ = — m{2n + 1). (29) 

To leading order in c , 

w = —i(n5 + (/)) — mcD, (30) 
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where <5 = 2n/S , 4> = (Ao^a + 7r)/<5o, Cj — 6 m /5o- So far, our results have been completely general and 
hold for arbitrary dimensions. It has been pointed out that the integrals for 60,6m can be expressed 
in terms of elliptical functions for D = 4 [13]. From eqs (26) to (28), it is easy to see that similar 
situation also occurs for D = 5 and 7. Therefore, we will only focus on these cases hereafter. 

Since no explicit expressions for the D = 4 case are given in Ref. [13], we think it is appropriate 
to present them here. The explicit form of qo(r) is given by 
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(31) 



There are four transition points. To leading order in c 1 , they are r = 0, r — — 2u\, and r = u\ ± ivi, 
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Here, E(m),E((p, m), K(m), F(ip, m),H(n, m) and U(n, ip, m) are the elliptical integrals: 
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Note that because of the pole at r — tq and branch cuts, there is ambiguity in eqs (37) and (38). To 
find the correct analytic expression, we compare the result with the numerical integration. After some 
trial and error, we find 
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where 9{x) is the step function. The term with a step function is introduced to compensate the 
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given by the residue of q^ 2 / A at r = r+. Similarly, 
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The magnitude of jump is 
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Note that here there is no discontinuity and we do not need to introduce any term with step function. 
Comparisons to the numerical results from Ref. [11] are given in Fig. 2 and 3. 
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Figure 2: Comparison between the analytic and numerical results for the real part of the highly 
damped QNM frequency u>(a). 
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Figure 3: Comparison between the analytic and numerical results for the level spacing 5. 
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In the extremal limit, a — ► /j,/2. As a result, A — > ^±v^,«i 
L'Hospital's rule, we find 
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The limit is smooth in constrast to the Reissner- Nordstrom case as pointed out in Ref. [13]. In the 
Schwarzschild limit a — > 0, even though the asymptotic QNMs are not continuous the level spacing 
does goes to the Schwarzschild result 5 = (D — 3)/(2ju) smoothly. 
For D = 5, 
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There are four transition points: r = ±.u\ ± ivi, with 
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Here, we have 
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The results are shown in Fig. 4 and 5. In contrst to the D — 4 case, there is a critical value 
a c Cl.678^ 1 / 2 . When a > a c , a) changes sign. The physical meaning is not clear at the moment. It 
might be related to the "algebraically special" frequencies [16]. 
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Figure 4: The analytic result for the real part of the highly damped QNM frequency o)(a) when D = 5. 
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Figure 5: Analytic result for level spacing 5 when D = 5. 

For D = 7, 

qo(r) =r 4 + a 2 r 2 + ^ia 2 r^ 2 . (55) 
There are six transition points: r = ±u\ ± iv\ and ±iw\. It is convenient to introduce the variable 
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The magnitude of jump here is again determined by the residue of q /A at r = r+. The results are 
shown in Fig. 6 and 7. There is no upper bound in a. Again, there is a critical value a c ~ 0.516/W 4 
and w changes sign when a > a c . 
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Figure 6: The analytic result for the real part of the highly damped QNM frequency u)(a) when D = 7. 



3 Conclusion 

In sum, we have calculated the asymptotic form of quasinormal frequencies to leading order for Kerr 
black holes in four, five and seven dimensions. All the releveant quantities can be expressed in terms of 
elliptical integrals. In particular, the results are consistent with the numerical ones in four dimensions. 
It would be interesting to generalize the method to other rotating black holes, especially Kerr-AdS 
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Figure 7: Analytic result for level spacing 5 when D = 7. 



black holes in five dimesions. It is pointed out in Ref. [17] that a large black hole in AdS background 
corresponds to a thermal state in CFT. Therfore, studying the decay of the scalar field may tell us 
something about the decay of a perturbation of the thermal state [18]. Extension to next to leading 
order is also desirable. One need to do systematic expansion in c _1 in a way similar to that used in 
Rcf. [7] . Such calculation may allow us to do a more comprehesive comparison between the analytic 
and numerical results. 
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